Revised value of the eighth-order electron g — 2 
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The contribution to tlie eightli-order anomalous magnetic moment (g — 2) of the electron from a 
set of diagrams without closed lepton loops is recalculated using a new FORTRAN code generated 
by an automatic code generator. Comparing the contributions of individual diagrams of old and 
new calculations, we found an inconsistency in the old treatment of infrared subtraction terms in 
two diagrams. Correcting this error leads to the revised value —1.9144 (35)(Q/7r)* for the eighth- 
order term. This theoretical change induces the shift of the inverse of the fine structure constant by 
-6.41180(73) X 10"^ 

PACS numbers: 13.40.Em,14.60.Cd,12.20.Ds,06.20. Jr 



The anomalous magnetic moment (.g—2) of the electron 
hasplayed a central role in testing the validity of QED 
P,[2. Recently, a Harvard group measured the electron 
g — 2 value Q using a Penning Trap with a cylindrical 
cavity [1]. Their result for ae = (.g-2)/2 3] 

ae = 1 159 652 180.85(0.76) x lO"^^ [o.66ppb], (1) 

has a 5.5 times smaller uncertainty than the best previous 
measurement 5j. 

To match the precision of this measurement the theory 
of g—2 must include up to the eighth-order contribution of 
the QED perturbation theory 0, i, 0, B B [lOl] as well as 



sented by self-energy-like diagrams related by the Ward- 
Takahashi identity. Taking the time reversal symmetry 
of QED into account, 518 vertex diagrams are amalga- 
mated into 47 self-energy-like diagrams shown in Fig. [1] 
Their integrands were carefully analyzed and checked by 
various means. However, no independent check of calcu- 
lation has been attempted until now. 



the hadronic [11|, [IJ] and weak contributions [13| . The 
tenth-order contribution of QED A^^^\a/'i:Y' might be 
relevant, but at present it is not known. As a temporary 
measure we adopt the bound A^^*^^ = (3.8) proposed 







y-fOr, y/I J\ ^77?i\ ^f1^. y^CT^ 


MOl 


M02 


M03 M04 M05 M06 M07 

^\ rCTy fCTx^ 






MIO Mil M12 MI3 MI4 

fC?^ /r?f>^i /fxj^ 


MI5 


M16 MI7 MIS MI9 M20 M21 

/tIC^x /f)C5\\ /rfr>v\ fjC^ 


M22 


M23 


M24 M25 M2fi M27 M28 

.77^ 


M29 


M30 


M31 M32 Mi3 M34 M35 

f^T^. /r^^ ff^^ /Y^\ 



in Ref. [ij to indicate a likely range of value taken by /7^^ frr:^\ 



A^j^°\ This will soon be replaced by an actual number 



which is being evaluated right now 

[3, liiH- Until 
then, the tenth-order term is the source of the largest 
theoretical uncertainty of the electron g — 2. The next 
largest uncertainty comes from the numerical integration 
of the eighth-order coefficient A^f^ [l3| • 

The purpose of this letter is to report the new value 

4.914 4 (35) (2) 



obtained by combining the information derived from the 
previous result [l^ and a new and independent evalua- 
tion of Af^ by means of FORTRAN codes generated by 
an automatic code generator "gencodeN" 15, ijj. 

A\ receives contributions from 891 Feynman dia- 
grams. 373 of them that have closed lepton loops had 
been calculated by more than two independent methods 
[lo| . The remaining 518 diagrams that have no closed 
lepton loop (called q-type) form one gauge invariant set 
(Group V) . In our formulation these diagrams are repre- 



FIG. 1: Eighth-order Group V diagrams. 47 self-energy-like 
diagrams of Moi - M47 represent 518 vertex diagrams. 

Technical progress in handling UV- and IR-divergences 
has enabled us to generate the A^th-order FORTRAN 
code easily and swiftly [H, Although "gencodeN" 
was developed primarily to handle the tenth-order term, 
we have applied it to fourth-, sixth- and eighth-order q- 
type diagrams as part of the debugging effort. With the 
help of "gencodeN" eighth-order codes are generated eas- 
ily. However, their numerical evaluation by VEGAS [l^ 
is quite nontrivial and requires huge computational re- 
source. Numerical work has thus far reached relative un- 
certainty of about 3 % . Although this is more than an 
order of magnitude less accurate than the uncertainty of 
the old calculation [l3|, it is good enough for checking 
algebra of the old calculation. 

Ultraviolet (UV) divergences of vertex and self-energy 
subdiagrams are removed by the if-operation 15|, Il9l . 
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TABLE I: Comparison of the numerical calculation of A/01-M24 of the eighth-order Group V diagrams. The second column 
shows the analytic expression of the difference of old and new calculations of the magnetic moment. The third column, value A, 
is obtained by plugging lower-order renormalization constants, such as AM4a, ALas into the expression in the second column. 
The fourth column, value B, lists the numerical values of AM°''* - AM"°". The fifth column is the difference A - B. If both 
numerical calculations are correct, A — B must vanish within the numerical uncertainty. In evaluating AM'^'^^ the double 
precision is used for the diagrams without a self-energy subdiagram, while the quadruple precision is used for the reminder. 



Diagram 


difference 


value A 


value B 


A- B 


Moi 








-0.0129(47) 


0.0129(47) 


M02 


2ALfifi Mr> 


-0.0066(3) 


0.0018(127) 


-0.0084(127) 


Mm 


ALRf-iMo 

'—^-'-'O J 0-^'-'- Z 


-0.1132(2) 


-0.1055(100) 


-0.0076(100) 


Mq4 


2fALK^i + ALr^-i^Mo 


0.3338(6) 


0.3515(221) 


—0.0177(221) 


Mo5 








0.0020(28) 


-0.0020(28) 


M06 








-0.0223(61) 


0.0223(61) 










—0 0102('40l 


0102f40l 


M08 


2(A5mi„AMA -1- ALinAMi„) 


-2.1809(7) 


-2.1773(163) 


-0.0011(163) 




lALRfyMo 


0.0805(2) 


0.0912(122) 


-0.0106(122) 


Mio 


2{ A&mAhAMA (^^ \ + ALRrioMo + ALa^AMa^.^ 


15.8899(49) 


15.8615(210) 


0.0284(216) 


Mil 


2ALfirfr,M9 


0.6948(3) 


0.6827(112) 


0.0121(112) 


M12 


2{ALeai + ALea3)M2 


1.2841(0) 


1.2875(74) 


-0.0034(74) 


Mi3 


2AL6hiM2 


-0.4202(4) 


-0.4238(48) 


0.0036(48) 


Mi4 


2AL695M2 


0.0892(3) 


0.0960(95) 


-0.0068(95) 


Mi5 


2AL691M2 


0.0889(3) 


0.0893(71) 


-0.0004(71) 


Mie 


2(A(5m4aAM4b(i.) + AL6C1M2 + ALisAMia) 


-2.6042(6) 


-2.6316(235) 


0.0274(235) 


Ml 7 


2(AL6el + AL6d4)M2 


-2.1183(5) 


-2.1010(189) 


-0.0173(189) 


Mis 


2{A(5m4i,AM4i,(i.) + ALisAM^b + {ALeti + AL6a2)M2} 


16.9690(39) 


17.1897(206) 


-0.2207(210) 


Mi9 








0.0002(3) 


-0.0002(3) 


M20 








0.0010(17) 


-0.0010(17) 


M21 








0.0003(3) 


-0.0003(3) 


M22 








-0.0090(25) 


0.0090(25) 


M23 


2AL6h2M2 


0.0495(3) 


0.0438(59) 


0.0057(59) 


M24 


2AL692M2 


0.0786(2) 


0.0945(61) 


-0.0158(61) 



120 . |21|. which is identical with the old approach. For 
diagrams containing self-energy subdiagrams, however, 
"gencodeN" treats UV-finite parts of self-energy subdi- 
agrams and IR divergences differently from the old ap- 
proach [1^ . 

Comparison of the new (still tentative) and old calcu- 
lations has revealed an inconsistency in the treatment of 
the infrared (IR) divergence in the latter, which is cor- 
rected in this letter. Thus we now have two independent 
evaluations of A\ . Of course, much more numerical 
work is required to reach the precision comparable to 
that of the old calculation. Fortunately, correction terms 
themselves can be evaluated easily and very precisely as 
are shown in ([4]) and 

Finite integrals AM°''^, z = 01, • • • , 47, from the previ- 
ous calculation are given in Ref. [l3|- AMf°'" are calcu- 



lated using the programs generated by "gencodeN" [15 , 
[l6l |. The numerical values corresponding to AMf'*^ — 
AM^™ are shown as value B in Tables |T] and |TT1 Since 
the diagrams without self-energy subdiagrams do not 
have IR divergence, AMf'^ and AM^""" should be iden- 
tical. This is confirmed within the numerical precision 
of AMf On the other hand, diagrams containing 
self-energy subdiagrams have IR divergence. The new 
treatment of their contributions produces results different 



from those of Ref. 0. The difference AA/f'^ - AMf"'' 
is listed symbolically in the second column of Tables U 
and [Til Their numerical values are calculated using the 
lower-order renormalization constants in Tabic IIIII and 
are shown as value A in Tables [J and [III The difference 
of value A and value B is listed in the fifth columns of 
Tables [J and [TTl If both calculations are free from error, 
value A and value B must agree with each other. 

Tables HI and HD show that "old" and "new" calculations 
are in good agreement for most diagrams. However, a 
large discrepancy —0.221 (21) is found for the diagram 
Mig. Though no detectable discrepancy is found for Afie, 
it has a structure similar to Afig and is somewhat simpler 
to analyze. Thus we examine here Afig instead of Afig. 

After an intense scrutiny of the programs of AM°q'^ 
and AAf "g*, our attention was focused on one of the IR 
subtraction terms of the finite term AM°g'^ 



211,123: 



AA/°e 



Id 



M- 



16 



j2 n I^-^^ie 

hclM2 - -JqcM2 - hsAMia 



ASniiaL 



4b(l' 



l2'ASm4aM2 



(3) 



where Miq is the bare amplitude, J^f Yisef ^sMi4 are 
the UV counter terms defined by the A'-operations 15|, 



3 





TABLE 11: Comparison of the numerical calculations of A/as- 


M47 of the eighth-order Group V dia^ 


rrams. 


Diagram 


difference 


value A 


value B 


A - B 


Mas 








-0.0031(20) 


0.0031(20) 


M26 


A5m6/(Ma. - Ma* [/]) 


2.5114(4) 


2.5369(95) 


-0.0255(95) 


M27 


2AL694Ma 


-0.0629(2) 


-0.0459(90) 


-0.0170(90) 


Mas 


2{A(5m6d(Ma. - Ma- [/]) + AL6c2Ma} 


-7.5329(6) 


-7.5310(189) 


-0.0020(189) 


Ma9 


2AL6e2Ma 


-0.2856(3) 


-0.2809(109) 


-0.0047(109) 


M30 


A5m6a(M2. - Ma* [/]) + 2AL662M2 


0.2768(7) 


0.2490(188) 


0.0278(188) 


M31 








0.0007(5) 


-0.0007(5) 


Maa 








-0.0024(10) 


0.0024(10) 


M33 








0.0001(3) 


-0.0001(3) 


M34 








-0.0010(13) 


0.0010(13) 


M35 








0.0001(13) 


-0.0001(13) 


M36 








-0.0027(22) 


0.0027(22) 


M37 








0.0004(5) 


-0.0004(5) 


Msg 


A5m6h(A'/a* - Ma* [/] ) 


-0.9088(3) 


-0.9112(40) 


0.0025(40) 


M39 








-0.0031(18) 


0.0031(18) 


M40 


2A5m6g{M2' - Ma. [/] ) 


3.8271(6) 


3.8326(71) 


-0.0055(71) 


M41 


A5m4a(AM4„(a.)) + AL4a;AM4a 


0.9809(3) 


0.9713(83) 


0.0096(83) 


M4a 


A5m6c(Ma. - Ma- [/]) + AL4! AM4a 


-7.0216(5) 


-7.0202(114) 


-0.0014(114) 




lAC fATl/T AC /TIT Tl/TfrlXI 

+Adm4a|AM4j,(a.) - Adma-(Ma. - Ma* 








M43 




0.4719(2) 


0.4703(42) 


0.0016(42) 


M44 


2AL6g3Ma 


-0.0751(2) 


-0.0499(69) 


-0.0253(69) 


M45 


A5m6e(Ma. - Ma* [/]) + ALecsMa 


-0.0515(6) 


-0.0498(90) 


-0.0017(90) 


M46 


A5m46AM4„(a.) + ALsesMa + ALi^AM4b 


-7.9336(22) 


-7.9232(86) 


-0.0104(89) 


M47 


ASmeb{M2- - Ma* [/]) + ALebsMa + AL4iAM4b 


10.5868(15) 


10.5864(102) 


0.0004(103) 



+A5m46{AM4i,(a.) - A5ma*(Ma. - Ma- [/])} 



lig . |20|, |21[, and the remainder are the IR subtraction 
terms. By a term-by-term comparison, we found finally 
that the IR subtraction term /4b(i.) was the culprit. 

Separation of an IR divergent part and a finite part of 
an integral is arbitrary. However, we must keep track of 
what is treated as the IR divergent part. In particular 
the IR subtraction term in AMi and one used to calcu- 
late the residual renormalization must be identical. All 
IR subtraction terms are summed up in the end, which 
gives a finite contribution as a part of the residual rcnor- 
mafization [U [H, HI] • What we found is that old FOR- 
TRAN codes of Iib{i') have different forms in AAfig and 



46(1*)- 



If we use /4b(i.) defined in Ref. [2l| as a part of 



AA^, 



4b(l-) 



we must add the correction term 



AM 



add 



-2 X 



{dz] 



G- 



4 J U2Y4 

X 22^2(1- Ai)3(l 

0.029 437 8 (98) 



(4) 



to AMffl^. The functions A„U,V in Eq. (g]) are defined 
in the I1237 limit of the diagram Afi^. For precise defini- 
tions of these functions see Refs. 



/Wig. tor prec 



21, 23|. The 



overall factor 2 comes from the time-reversed diagram. 
The value ^ is smaller than the uncertainty of value B 
for Afig. Thus it is undetectable by direct comparison of 
values A and B until precision of AMfg™ is improved. 
Analyzing the difference of Afjg^ and Af"g in the same 



manner, we found that the correction term is not small 
for Afig: 



AAf- 



add 



= -2x ^ J{dz)G{l-K5) 

^ I U2vi ^2^2(1-^1)^(1-^2) 
= -0.215 542 (19), 



(5) 



where all Ai,U,V are defined in the I1237 limit of Mis- 
Their explicit forms are different from those of Afig. The 
function Sm4a(b) [fo] in -^16(18) ^® related to the UV-finite 
part A6m4a{b) of the mass-renormalization constant. If 
we add AA//f|'^ to AMf^'^, value B of Afig becomes 
16.974 (21) and the difference between values A and B 
is reduced to —0.005 (21), which is consistent with zero 
within the precision of numerical calculation. 

We should like to emphasize that the development of 
automatic code generator 1^ 1^ was crucial in discov- 
ering the existence of extra IR subtraction terms in A/ig 
and Mig. Details of our investigation will be reported 
elsewhere Adding the terms Eq. ( HI) a nd Eq. ^ to 
the "old" calculation Eq. (58) of Ref. [l3|, we find the 
entire contribution of Group V: 



Af^(GroupV) = -2.179 16 (343), 



(6) 



which is in good agreement with the still tentative value 
obtained by the code generated by "gencodeN" : 

Af)s™^(GroupV) = -2.205 (54). (7) 
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TABLE III: Finite renormalization constants used in Ta- 
bles U] and ini The validity of the sixth-order renormaliza- 
tion constants are checked by comparing the sum Xlbd = 
^Lfixi + ^ABqx + 2ASmex, x — a, ■ ■ ■ h to the previous 
Xlbd values listed in Ref. \Wk - 





539589f67~l 


A_/v6bl 


— 1 479347('244~) 




—0.167211(81) 


ALk^o 


0.582944(106) 




1 489038('142"l 




—0 016344f73~) 


A Lfi^i 


—0 21Q31 1 048^ 




—0 740890^373") 




OTlBWflSSl 




—0 285566(252^ 




—0 551410f236l 




—0 141327(380^ 


'— »-^Dal 


833454f402l 




088899(251^ 


ALk^o 


-0.090653(141) 




0.078625(184) 




—0 499683 ('4071 




—0 075127(176^ 


A fiR^A 


— 1 3774.50^287^ 


A T/K /I 


—0 062906(155") 


A ^jK^S 


69483.5 ^227^ 


A T,R„r 


089234(288^ 


A/!/f;fi 


—0 006638f212l 


ALqIiI 


—0 420233(330) 


A /!/K f o 


080534('139l 


ALq}i'2 


049517(284^ 


AZ/6/3 








ASm^a 


-0.15309(34) 


A5mgi, 


1.83775(25) 


ASniec 


-3.05039(22) 


ASmfid 


-1.90114(15) 


ASniee 


0.11210(25) 


A5m(,f 


1.25572(19) 


ASmsg 


0.95677(13) 


ASmeh 


-0.45439(11) 


ALic 


0.003387(16) 


ALi^ 


-0.481834(54) 


ALis 


0.407633(20) 


ALu 


0.124796(67) 


ASniia 


-0.301485(61) 


A5m4i, 


2.20777(44) 


AM^a 


0.218359(39) 


AM46 


-0.187526(39) 


AM4,(i.) 


3.61946(83) 


AM4a(2.) 


-3.60244(67) 


AM46(i.) 


4.25054(23) 


AM46(2.) 


1.64475(10) 


AMa 


0.5 


AM2. 


1 


AM2' [I] 


-1 


A5m2» 


-0.75 



The revised contribution ([6]) shifts the total eighth- 
order term ^4^*'' to the one given in Eq. As a con- 
sequence, the theoretical prediction of ae is moved by 
-5.421 775 (62) x IQ-^^^ yielding 



ae(Rb) 
ae(Cs) 



1 159 652 182.78 (7.72)(0.11)(0.26) x 10"^^ 
1 159 652 172.98 (9.33)(0.11)(0.26) x 10"^^ 

(8) 



where 7.72 and 9.33 come from the uncertainties of the 
input values of the fine structure constant 

a-\Rh06)= 137.035 998 84 (91) [6.7ppb], (9) 
a-i(Cs06) = 137.036 000 00 (110) [8.0ppb], (10) 

determined by the Rubidium atom [13] and Cesium atom 
25I 2^ experiments, respectively. The uncertainty 0.11 
of Eq. ([8]) comes from the eighth-order calculation and 
0.26 is an estimated uncertainty of the tenth-order term. 

Because of high precision of the experiment ([T]) the fine 
structure constant a determined from the theory and 
the measurement is sensitive to the revision of theory. 
The inverse fine structure constant a~^(ae) moves by 



-6.411 80 (73) X 10-'' from the previous value in Ref. 
The revised Q;~-^(ae07) is about 4.7 ppb (or about 7 s. d.) 
smaller than Q;~^(ae06), but is still in good agreement 
with a-i(Rb06) of Eq. (P and a-i(Cs06) of Eq. (fTOl) . 
whose uncertainties are about 7 ppb. 
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